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Abstract. The convex hull property is the natural generahzation of maxi- 
mum principles from scalar to vector valued functions. Maximum principles 
for finite element approximations are often crucial for the preservation of qual- 
itative properties of the respective physical model. In this work we develop a 
convex hull property for Pi conforming finite elements on simplicial non-obtuse 
meshes. The proof does not resort to linear structures of partial differential 
equations but directly addresses properties of the minimiser of a convex en- 
ergy functional. Therefore, the result holds for very general nonlinear partial 
differential equations including e.g. the p-Laplacian and the mean curvature 
problem. In the case of scalar equations the introduce techniques can be used 
to prove standard discrete maximum principles for nonlinear problems. We 
conclude by proving a strong discrete convex hull property on strictly acute 
triangulations. 



1. Introduction 

Let the bounded polyhedral domain cz M", n 5= 2, be partitioned by a con- 
forming simplicial triangulation T. We denote by V(T)™, m e N, the finite element 
space of m-dimensional vector valued, continuous, and piece wise affine functions 
over T and by Vo(T)™ c V(T)™ the subspace of functions with zero boundary 
values. For G e V(T)'" we consider minimising problems of the form 

(1.1) J{V) := [ F{x,\\7V{x)\)dx^mm inG + Vo(r)". 

Jo 



We suppose that (1.1) admits a unique minimiser and that the mapping fl x Mj> 3 
(x, y) t-^ F{x, y) e M is monotone in its second argument; for the precise assump- 
tions see Section [221 



The following theorem is the main result of this work. 

Theorem 1.1 (Convex hull property). Let the triangulation T he non-obtuse, i.e., 
the angle between any two sides of any simplex T e T is less or equal 7r/2. 
Then the unique minimiser U e G + Vo(T)"' of ( |1.1[ ) satisfies 

(CHP) U{n) c convhuU {U{dn)) = convhull {G{dn)). 

This is the discrete analogue of the so called convex hull property of vector 
valued minimisers; see |DLM981 IBF02) . Note that the convex hull property is the 
generalization of the maximum principle to the vector valued case. Indeed, if m = 1, 
then the convex hull property reads as 

U{n) = [mmU{dn),meixU{dn)]. 
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The problem above is presented as a minimising problem. If the energy J is 



Frechet differentiable, then (1.1) becomes equivalent to finding a solution U e 



G + Yq (T) of the Euler-Lagrange system 

(1.2) (DJ{U), H)= I a{x,\S/U{x)\)S/U:\7Hdx = for all e Vq (T ) , 
Jn 

where a{x,t) = Ft{x,t)/t. Consequently, Theorem |1.1| applies to a large class of 
non-linear partial differential equations. 

The research on discrete maximum principles can be traced back to the early pa- 
per of Ciarlet and Raviart |CR73| in the nineteen seventies. Provided a non-obtuse 
triangulation, they established a discrete maximum principle for scalar valued lin- 
ear second order elliptic problems on a two dimensional domain. This result was 
generalized to three dimensional domains in |KQ95, . For the Laplacian, the above 
mesh condition can be replaced by a weaker, so called Delaunay condition. This 
condition only requires that the sum of any pair of angles opposite a common side 
is less or equal tt; see [SF731 lLet92j . Draganescu, Dupont, and Scott [DDS05j ana- 
lyzed the failure of the discrete maximum principle when the Delaunay condition is 
violated. However, according to |RS821lKKN01| a discrete maximum principle may 
hold even if both angles in such a pair are greater than tt/2 introducing additional 
restrictions involving a larger neighborhood. Some of these concepts are generalized 
to anisotropic diffusion problems in jLHlOi IHuallj . 

Recently, discrete maximum principles are proved for quasi-linear partial differ- 
ential equations |KK05ai IWZ12j , which means in the case of ( |1.2[) that < A < 
a{x, t) ^ A < GO for all X, t e O X M. Jiingel and Unterreitner [JU05' proved a 
discrete maximum principle for partial differential equations with non-linear lowest 
order term. 

All aforementioned approaches are based on the ellipticity and continuity of 
the problems in order to allow for a Hilbert space setting. In particular, beside 
[WZ12[|JU05| they rely on the inversion of so called Af-Matrices. As a consequence 
the results are basically restricted to energy functionals with quadratic growth. 

In contrast, our approach is based on a projection into a suitable convex set 
and direct properties of the energy functional. Therefore, it does not rely on linear 
structures and directly applies to general convex non-linear, vector valued problems. 
As a drawback, we need point-wise (resp. element-wise) properties of the finite 
element functions and therefore cannot weaken the mesh restrictions to a Delaunay 
type condition. 

The plan of the paper is as follows. In Section [2] we introduce the finite element 



framework and provide precise assumptions on the energy functional from (1.1 ). In 
Section [3] we prove the discrete convex hull property. Theorem |1.1[ We close the 
article in Section [4] discussing several applications and extensions of the introduced 
theory. In particular, we prove a maximum principle for non-linear scalar partial 
differential equations with non-positive right hand side and show that the p-Laplace, 
as well as the mean curvature problem fit into the presented theory. Moreover, we 
show a discrete maximum principle for non-linear problems with lower order terms 
using mass-lumping. Finally, we define and verify a strong convex hull property, 
which generalizes the strong maximum principle of scalar valued problems. 

2. Problem Setting 

In this section we shall introduce the finite element setting and formulate the 
precise assumptions on the non- linear energy. For the sake of presentation we denote 
the Euclidean norms in M", n 5= 2, M™, m e N, as well as the absolute value in M 
by |-|. Furthermore, for A = [A,j] ,B = [B,.j] .-i „ e W^'""'' we define the 
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inner product A : B := AijBij and denote also by |A| the corresponding 

J=1,...,TT1 

Frobenius norm of A. 

2.1. Finite Element Framework. Let T be a conforming partition of the poly- 
hedral domain Q, a M" into closed n-simplexes. To be more precise, we have 

^ = [j{T\T eT) 

and the intersection of two different elements in T is either empty or a complete 
fc-dimensional sub-simplex, ^ fc < n. Let Pi(T) be the space of afhne linear 
functions on T e T and define the space of first order Lagrange finite elements by 

(2.1) V(r) := {V e C(f7) : Vjr e Pi(T) for all T e T} . 

This space is spanned by the so called nodal Lagrange basis. To be more precise, 
let Af be the set of vertices of all elements in T and let A/q := A/" n be the subset 
of vertices inside il. For z e Af the corresponding Lagrange hat function (/f^ is 
uniquely defined by (/i^ (y) = Szy for all y e J\f, where we make use of the Kronecker 
symbol, i.e., S^y = 1 whenever z = y and S^y = else. We denote the set of element 
sides of T by 5. For two sides of S'l, 5*2 G 5 of a simplex T e T we define the angle 
between 5*1 and S2 by Z(S'i, 5*2) := tt — Z(ni, 712), where Z(ni, 712) is the angle in 
between the normals rii and n.2 pointing inside T from Si and ^2, respectively. 

In the following we introduce non-obtuse meshes, which play a crucial role in the 
analysis of discrete maximum principles; see e.g. |CR73[ |KQ95[ [Dl)S05j . 

Definition 2.1 (non-obtuseness) . A simplex T e T is called non-obtuse if the 
angles between any two of its sides are less or equal ^ . 

The conforming triangulation T of Q is called non-obtuse if all simplexes T e T 
are non-obtuse. 

We remark that in higher dimensions n 5= 3 the construction and refinement of 
non-obtuse partitions is a non trivial task. The interested reader is e.g. referred to 
|KQ95llKKn5bllKKllj . 

One can characterize non-obtuse meshes by means of the nodal Lagrange basis 
functions. Indeed, observing that the gradient of a Lagrange basis function 
z e Af ri T on a simplex T e T points in the direction of the normal of the side 
opposite z, one can easily show that the following condition on the basis functions 
is equivalent to the geometric mesh property in Definition |2.1[ compare also with 
[CR73l|KQ95l . 

Proposition 2.2. An n-simplex T e T is non-obtuse if and only if 

V(j)z\T • '^4'y\T ^ for all z,y e T n Af with z ¥= y. 

A direct consequence is the following corollary. 

Corollary 2.3. Let T be a conforming triangulation offl. Then T is non-obtuse 
if and only if 

■ V(/>y ^0 a. e. in VL for all z,y e Af with z ¥= y. 

The subspace of functions in V(T) vanishing at the boundary dil is denoted by 

Vo(r) := {V e V(r) : V\en = 0} = span {0, : z e A^o} . 

For m e N, the m-dimensional vector valued finite element spaces corresponding to 
V(T) and Vo(T) are given by 

m m 

v(r)" := X v(r) and Vo(r)" := X Mr), 

i = l i=l 
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respectively. Clearly, for G e V(T)™, the affine subspace of functions that coincide 
with G ononis G + Vo(r)™. For uj Cl we define V{u}) := {V{x) : xeu} czW 
as the set of function values oi V on w. 

2.2. Energy Minimisation. In this section we formulate standard conditions on 



the energy in order to guarantee the unique solvability of (1.1). To this end we 
first ensure that JiV) is well defined for all V e V(T)'" assuming that i^: fixK^. 
M defines a Nemyckii operator T : V(T)'" L^(ri) setting 

{TV){x) ^ F{x,\VV{x)\), xen. 

This can be achieved by the following conditions. 

Assumption 2.4 (Nemyckii Operator). We assume that F: il x M. 

(1) is a Caratheodory function, i.e., F(-,t) is measurable in f2 for each t e 
and F(x, •) is continuous on for almost every x e Q; 

(2) satisfies the growth condition 

\F{x,t)\^a{x) + -f{t), {x,t)eVlxM.^^ 

for some a 6 L^{n) and continuous 7 : — > M^. 

Proposition 2.5. Suppose that F: ft x Mj. R satisfies Assumptions 2.4 Then 
the energy functional J : V(7~)™ — >■ -L^(il) defined in (1.1 1 is continuous. 



Proof Let {VfejfeeN c V(r)", V e V(r)™ such that Vk ^ V in V(r)" as 
k ^ CO. The space V(T)™ is finite dimensional and hence all norms on V(T)'" 
are equivalent. Therefore, we can conclude that the sequence {V^lfegN is uni- 
formly bounded in L'^^fl)"^. Assumption 2!^ |2) thus implies for some C > 

1) it follows that 



2.^ 



that Vfe(-))| ^ a(-) + C e L^{n) and from Assumption 

Vk{-)) F{-, V{-)) almost everywhere in fi. Hence the claim is a consequence 
of Lebesgue's dominated convergence theorem. □ □ 

The following conditions are linked to the existence of a unique minimiser of J^. 

Assumption 2.6. We assume that F : fl x Mj. — > M is strictly convex in its second 
argument, i.e., it holds for all 9 e (0, 1) and s,t^O with s =^ t, that 

(2.2a) F{-,es + {l-e)t) <eF{-,s) + {l-9)F{-,t) a.e.mQ. 

Moreover, assume that F satisfies the coercivity condition 

(2.2b) F{-, t) ^ g{t) for all t ^ Q a.e. in 

with some continuous monotone function g : — > M such that g(t) — > <X as i —> co. 



Now we are prepared to prove the existence of a unique solution to ( 1.1 1 



The 



Proposition 2.7. Let G e V(T)™ and suppose Assumptions 2.4 and 2.6 
there exists a unique U e G + Vo(T)™ such that 

J{U) = min {J{V) : y e G + Vo(r)"} . 

Proof. The proof of this assertion is fairly standard. However, for the sake of a self- 



consistent presentation we sketch the proof here. Thanks to ( 2.2b ) it follows that the 
energy functional is bounded from below and that an infimal sequence {VfelfceN c 
G-(-Vo(T)™ is norm-wise uniformly bounded by some constant C > 0. Since V(T)'" 
is of finite dimension, { V^jfegN is pre-compact in V(T)'". Consequently, there exists 
a converging subsequence Vk,, —* U for some f/ e G + Vo(T)'" as ^ ^ co. By the 



continuity of J , Proposition 2.5 it follows that J7 is a minimiser of J . Finally, the 
strict convexity (2.2a) implies the uniqueness of the minimiser. □ □ 
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Figure 1. The projection Uk 



The following assumption is crucial in order to prove the discrete maximum 
principle (Theorem |l.l[ ). 

Assumption 2.8 (Monotonicity). We assume that the mapping F{x, •) : — > M 
is monotone increasing for almost every x e Q. 

Note that V(T) c L^(f2) is finite dimensional. Therefore, in Assumptions 2.4 
and |2.6| we do not need to resort to common growth or coercivity conditions like 
7(t) = c|i|^ or g{t) = c|i|'' for some c > and s > 1. As a consequence, the 
presented theory can be applied e.g. to the mean curvature problem; compare also 
with Remark 14.31 below. 



3. Proof of the Discrete Convex Hull Property, Theorem |1.1| 

In the following let K a denote a convex, closed subset. We denote by 
Uk ■ R™ — » K the orthogonal projection with respect to the Euclidean inner 
product of . In other words 

(3.1) IIkx := argmin^^^ \x - y\. 

The next lemma is an immediate consequence of the convexity of the set K. 

Lemma 3.1. Let K cz M™ be a closed and convex. Then for all x e M™ and z e K 

we have 

{x - Ilif x) • (z - Wkx) ^ 0. 

Proof, li x e K, then x = Tlj^x and the estimate is obvious. Let us assume that 
X ^ K. We denote by the hyper-plane which touches K in Hxx and is orthogonal 
to a; — Hkx. Since K is convex, separates K from x. This implies that the angle 
between {Hkx — x) and (JIkx — z) is greater or equal 7r/2; compare with Figure [ij 
This proves the assertion. □ □ 

We define a projection operator Vk ■ V(r)" ^ V(r)" with ('PkV')(17) cz K hy 
setting 

(3.2) {VkV){z) ■.= UkV{z), zeAf. 

We next list some elementary properties of this projection operator. In particular, 
since K is convex and V(T)™ contains only piecewise linear functions over T, we 
have 

(3.3a) 



VkV -.n^K for all y e W{TY 



6 



L. DIENING, C. KREUZER, AND S. SCHWARZACHER 



If G{d^l) c K, for some G e V(T)™, it follows that the boundary values are 
preserved, i.e., 

(3.3b) VkV eG + Vo(r)™ for al\VeG + Vo(r)"\ 

The following property of Vk is the key estimate in the proof of the convex hull 
property. 

Lemma 3.2. Let T be a non- obtuse conforming triangulation ofD, and let K c M™ 
be a closed and convex set. Then for V e V(T)™ we have 

VV : VVkV ^ \VVkV\^ a. e. in VL, 

and 

\VV\ ^ |V-Pi^y| a. e. in Q.. 

Proof. We prove the estimates on a fixed n-simplex T e T. Let zq, . . . , denote the 
vertices of T and let 00, . . . , 0„ be the corresponding local Lagrange basis functions. 
Then 

n 

VlT = 2 V(Z,)0.|T 

1=0 

and 

n n 

Wit : V(PkV)|t = 21] (Viz,) ■ IlKV{zk)) (V(/),;|t • V0,|t) 

1=0 fe=0 

Using V0O + V(/ii + • • • + V(/)„ = VI = on T we conclude 

n 

(3.4) VFiT : y{VKV)\T = E E ^(^*) ' i^KVizk) - n;^y(z,)) V0,|T • V</)fc|T. 

i=0 k+i 

Lemma |3.1| yields 

{Viz,) - UKViz,)) ■ (UKVizk) - TlKViz,)) 0, 

and thus 

V{z,) ■ {UKV{zk) - UkV{z,)) < UkV{z,) ■ {UKV{zk) - UKV{zi)). 

Since T is non-obtuse we have by Proposition |2.2| that V4)i\T ■ '^4'k\T ^ for all 
i ¥= k. Consequently, we arrive at 

n 

VFiT : V(7'a'V)it 5= J] J] {IIkV{z,) ■ {nKV{zk) - IiKV{z,)) (V<^,it • ^^^kir) 

n n 

= 22 (nA'V(z,) • UKV{zk)) (V(/.,|T • V0fe|T) 
i=0 fc=0 

= V(PkV)|t : V(7'kV)|t = |V(PkV)|t|'. 

This proves the first claim. 

Thanks to the Cauchy-Schwarz inequality, the second estimate is an immediate 
consequence of the first one. □ □ 

Lemma 3.3. Let T be a non-obtuse conforming triangulation of and let F : 
nxR^ satisfy As sumptions \2^\2^ and\t^ For G e V(r)™, meN, let 

K cz M'" be a convex and closed set such that G{dQ) c K. Then the uniguely 
defined solution U e G + Vo(T) '" of satisfies 

U{n) c K. 
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Proof. It follows from Proposition 2.7 that problem (1.1 1 has a unique solution 
U e G + Vo(T)™. Since K is a closed convex set we conclude from Assumption 2.8 
and Lemma 13.21 that 

J{U) > JiVKU). 



It follows by ([3^ that VkU e G + Vo(r) 
of J in G + Vo(r)'" we conclude VkU = 



consequence of (3.3b) 



™. Since U is the unique minimiser 
U. The desired assertion is then a 
□ □ 



From this result the discrete maximum principle, Theorem |1.1[ follows as a 
corollary. 



Theorem The convex set K 
of Lemma [3.31 



convhuU(G((?rj)) 
□ 



satisfies all conditions 
□ 



Remark 3.4. We emphasize, that the proof of the convex hull property (Theo- 
rem \l.l\ does not directly resort to Assumpti 



lions 



2-4 and 2.6. which are only use 



in Proposition \2. 7| to ensure the existence of a solution. In fact, a unique min- 
imiser U e G + Vo(T) of ( 1.1 1 satisfies the convex hull property (Theorem \l.l\ l if 
Assumption 2.8 holds. The uniqueness of the minimiser is usually a consequence 
of the strict convexity of the energy. 



4. Applications and Extensions 

In this section we apply the above techniques to some particular exemplary 
problems and extend Theorem |1.1| in several ways. 

4.1. The p-Laplace Problem. This section is concerned with the so-called p- 
Laplace operator. We verify that p-harmonic functions fit into the framework of 
Section [2] and then prove a discrete maximum principle for scalar solutions of a 
p-Laplace problem with non-positive right-hand side. To this end let T be a con- 
forming triangulation of the n-dimensional domain D, with polyhedral boundary, 
G 6 V(T)™, TO e N, and p,qe (1, oo) with i -f 1 = 1. It is weU known, that 

(4.1) I \\/Uf^^\/U : \/V dx = in O 

uniquely determines a function U e G + Vo(T)™. The equivalent minimising 
problem reads as 

J{V) := \ -IVVf dx ^ min inG + VofT). 
Jn P 

It can easily be verified that the operator 

F{x,t) := -\tf, t^O 
P 

satisfies Assumptions |2.4[ |2.6[ and |2.8| Consequently, if T is non-obtuse, all condi- 
tions of Theorem 1 1.1 1 are satisfied. 

For TO = 1 and / e L''(fl) we investigate the non- homogeneous p-Laplace problem 

( |Vf7|^"^V[/- VV^da;= | fVdx in Q 

(4.2) Jo Jn 

U = G on dfl. 

The unique solution U e G + Vo(T) is the minimiser of the energy 

Jf{V):= \ -\VVf-fVdx, VeG + Yo{T). 
Jo P 
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Since the integrand depends not only on the gradient of the function but also on 
the function value itself, we can not directly apply Theorem However, by 

modifying the arguments we can prove the following result. 



Let the conforming triangulation T of Q be non-obtuse and let 



Theorem 4.1. 

f ^0 mQ. 

Then the unique solution U e G + Yo{T) of (4.2 1 satisfies 

(DMP) ma.xU{n) ^maxU{dn). 

Proof. We define the closed convex set 

K := {-00, max U{dn)] 3 conv hull {U{dn)) . 

Hence the projection operator Vk ■ V(T) V(T), defined in ( |3.2| ) satisfies (3.31 
i.e., we have that VkU e G + Vo(T). Clearly the mapping {x, t) i— > F{x, t) := 



HP, 



{x,t) e Q, X Kj., satisfies all conditions in Section 2.2 In particular, it satisfies 
Assumption |2.8[ i.e., it is monotone increasing in its second argument. Therefore, 



we obtain by Lemma 3.2 that 

f -IVVkUI" dx f -iVUl^dx. 
Jn P Jn P 

Moreover, by means of ( |3.1[ ) and the particular choice of the convex set K, we have 

VkU{x) = min{t/(a;),maxC/((5r2)} U{x), for a\\ x e Vl. 
Hence the assumption / ^ in implies 



(4.3) 



r fVxUdx^ [ fUdx. 
Jn Jn 



Combining these estimates yields 

JiVKU) < J{U). 



Since U is the unique minimiser we arrive with (3.3) at U{fl) = (VKU){fl) cz K. 
This proves the assertion. □ □ 

Remark 4.2 (Orlicz functions). The p-Laplace operator often serves as a prototype 
for a much larger class of non-linear operators: Indeed, if we have a continuous, 
strictly convex and monotone function : IRj> — > M^, then the non-linear function 
F{x,t) := ip(t) satisfies Assumptions 2.^. 2.6. and 2.8 Consequently, the discrete 



convex hull property (Theorem 1.1) holds also for this larger class of problems. 

In order to prove Theorem 4-1 for more general nonlinearities we need to ensure 
the solvability of the problem. This can e.g. be achieved either by restrictions on the 
right-hand side f or by assuming the additional growth condition limj^oo '4'it)/t = go 
for the non-linearity. We illustrate the first case in Remark \4-3\ below by means 
of the mean curvature problem, i.e., il'{t) = Vl + t^- The latter case includes 
for example all N-functions. For more information on N-functions and on finite 
element approximations for this kind of problems we refer to |DR07llDK08|,lBDK12j . 

Remark 4.3 (Mean curvature problem). We emphasize that also the discrete mean 
curvature problem fits into the setting of Section \2.S\ To see this, we consider the 
corresponding energy 



JiY) 



I 

Jn 



i + |vyrda;, veYiry 



Obviously, the mapping Q x IRj> 3 (x, t) i— > F(x, t) := \/l + t? has linear growth and 



is monotone in t. Therefore, it satisfies Assumptions 2.4 and 2.8 Moreover, it is 
coercive and thanks to ^ Vl + = ^^^^^2 > f^i^ t ^ 0, it is strictly convex. 
This verifies Assumption 2.6 and consequently, minimisers U e G + Vo(T)™ of the 
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energy functional J in G + Wo{TY'^ , G e V(T), satisfy the convex hull property 
Theorem li.il 

Consider the scalar case m = 1 and let f e L^(fl), f ^ in fl. Let G e V(T), 
according to |FV03j the energy 

(4.4) Jf{V) := f Jl+\WV\^-fVdx, VeG + Yo{T), 

Jn 

is coercive if and only if there exists e > 0, such that 

f fVdx^{l-e) r |Vy|da; for all V e YoiT). 
Jo Jn 

Under this condition Theorem \4-.l\ directly applies to minimisers of the energy de- 
fined in (4.4). 

4.2. Problems Involving Lower Order Terms. In this section we investigate 
non-hncar problems involving lower order terms in form of lumped masses. Let T 
be a conforming triangulation of f2 and G e V(T)™, m e N. For 1 < p, q < X) let 
the energy 

(4.5) J{V) := i I Vy r + i I Vrda; Ve V(r)"^. 

serve as our prototype. Obviously, there exists a unique minimiser of in G + 

Vo(r)'". 

In order to prove a convex hull property (Theorem in the spirit of Section [sj 
we would need an analogue of Lemma |3.2| for function values. Unfortunately, in 
general we do not have |V| 5= I'PxV'l point-wise in O. One can construct examples 
where the discrete maximum principle is violated even in one dimension; see e.g. 
[BKKOSj . However, for p = g = 2 it is well known that maximum principles hold for 
sufficiently small mesh-sizes if the mesh is strictly acute; see e.g. |CR73l [BKK08| . 

This additional assumptions can be avoided by mass-lumping, i.e., instead of a 
minimiser of ( |4.5[ ) we look for a minimiser of the modified energy 

(4.6) MV) := f i|vyrd.+ 2 Vwi', 

which is coercive and strictly convex. Hence there exists a unique U e G + Yq{T)"^ 
minimising J^i in G + Vo(T)™. This minimiser satisfies the following convex hull 
property. 

Theorem 4.4. Let the conforming triangulation T of Q be non-obtuse and let 
U e G + Vo(T)™ be the unique minimiser of the energy Ji in G + Vo(T)™. 
Then we have the max hull property 

U{n) cz convhun(Lr((517) u {0}). 

Proof. Thanks to K := convhull [G{dft) u {0}) we can choose z = in Lemma 
and observe that 

^ \Ukx\^ - X ■ Hkx ^ \Q.Kx\^ - \x\\X\kx\ for all x e M'". 

Consequently we have 

\IIkU{z)\ ^ \U{z)\ for all zeAf. 
Combining this with Lemma |3.2[ we obtain 

MVkU) ^ J,{U). 

Recalling that the minimiser U is unique in G + Vo(T)™ we have U = VkU . 
Hence we can conclude with (3.3a) that [/(fi) = {VkU){^1) c K. This proves the 
assertion. □ □ 



3.1 
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Remark 4.5. In the scalar case m = 1 the minimiser U e G + Wq{T) of the energy 

with / < satisfies the discrete maximum principle 

maxC/(f7) ^ niax{0, maxU{dfl)} . 



This claim is an easy consequence of the above arguments combined with (4.3). 



4.3. Strong Convex Hull Property. In some situations we have a strong version 
of the convex hull property. For scalar functions this is well known as the strong 
maximum principle: if a subharmonic function attains a global maximum in the 
interior of the domain fi, then the functions must be constant (on the connected 
component). 

We will introduce in this section a discrete strong convex hull property. To this 
end we need that the conforming triangulation T of is acute: A simplex T e T 
is called acute if the angles between any two of its sides are less than ^ . Similar to 
Definition |2.1[ we can formulate this definition also in terms of the Lagrange basis. 



Definition 4.6 (acuteness). Let T be a conforming triangulation of CI. We call an 
n-simplex T e T acute if 

V(j)z\T • '^4'y\T < for all z,y e T r\ Af with z ¥= y. 

The conforming triangulation T of CI is called acute if all simplexes T e T are 
acute. 

Let K c M™ be a non-empty, convex, and closed set. Then we call x e K an 
extreme point if it is no element of a line segment of two different points of K. By 
exii{K) we denote the set of extreme points of K. It is well known that if K is 
additionally bounded, we have conv hull(extr(i^)) = convhull(A'). 

We assume the following additional condition on our mesh. 

Assumption 4.7. We assume thatT is a conforming triangulation of CI such that 
for every T e T , we have T n Ao . 

In other words, every simplex of the triangulation T has at least one interior 
vertex. As a consequence, we have that every vertex of the boundary is in the 
support of one of the interior Lagrange basis functions. Assumption |4.7| is not a 
severe restriction. 

We can now state the discrete strong convex hull property for discrete p-harmonic 



functions, i.e., solutions of (4.1) 



Theorem 4.8. Let the acute conforming triangulation T of CI satisfy Assumption 



.1 Let U e G + Vo(T)™, be discretely p-harmonic with boundary values G 6 



V(T)'", i.e., U is a .solution o/ (4.1). 



IfU{zQ) e extr(conv hull(C/(il)) for some Zq e Cl, then U is constant. 

Before we get to the proof of Theorem |4.8| we need a local version of this result. 



Lemma 4.9. Let T be acute and let U be a solution of (4.1). Let zq e Ao and 

let C := (Af\{zQ}) n suppc/)^;,, i.e., C is the set of neighbors of zq. Then U{zq) g 
extr(convhull([/(C))) implies that U is constant on supp^z^. 

Proof. Let us assume that U is non-constant on supp (jjzo ■ We have to prove that 
U{za) i extr(convhuU(C/(C))). 
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Since zq e Aq, we have for every unit vector Bj e M™, j = l,...,m, that 
ejcpzo 6 Vo(T)™. Therefore, it follows from K^ that 



r \S/Uf ^V[/ : V(ej0^J dx = for all j = 1, . . . ,m. 
Jo 



Thanks to the representation U = XlyeAA ^iy)^yy thus get 



or equivalently 

U{zo) f |Vi7r^'V(^,„ • V0,„ = VUiy) f |Vt/r'(-V^y • V,^,„) dx . 
^ , yec h ^ . 

By assumption (7 is non-constant on supp^z^, hence /3o > 0. On the other hand, 
since T is acute, we have /Sy 5= for all ye C. Let T e T, T c supp such that 
f/|r #0. li y e Af n T, y Zq, then the acuteness of T implies /3y > 0. Thus at 
least two (more precisely: at least n) of the /?y with y e C satisfy (3y > 0. 

The identity (jj^g + Yjy^c't'v = 1 on supp^i^o implies /3o = Y^y^c Py Therefore, 
defining A^^ := Py/Pzo yields U{zo) = J^yec ^v^(y) "^i*^ '^v > for y e C and 
l^yeC ^ 1- Consequently, at least two of the Aj^ with z e C satisfy Ay > 
and hence U{zo) cannot be an extreme point of conv hull {J7(j/) : y e C}. In other 
words, we have U{zo) ^ extr(convhull(L''(C))). □ □ 

We can now prove Theorem |4.8[ 



Theorem \4.8\ Since U is piecewise linear and thanks to the properties of extremal 
points, it follows that there exists zq e Ao with U{zq) = U{zo). Therefore, we may 
assume w.l.o.g. that zg e Aq. 

Suppose U{zo) e extr(convhull(C/(ri)) and let C := (A/\{zo}) n svipp(j)zo be 
the neighbors of zq- Due to the convex hull property of U (see Theorem 
we have that conv hull(t/(C)) c: convhull([/(c/ri)). Since U{zq) is an extreme 
point of convhull(L''(c/f2)) it must also be an extreme point of extr(convhull((7(C)). 
Therefore, by Lemma |4.9[ U is locally constant on supp^zo- 

A repetition of this argument shows that U is constant on its connected compo- 
nent, which is 51 itself. For this argument we need, that every simplex contains at 
least one vertex of Aq. □ □ 



Remark 4.10. We emphasize that (4.1) serves as a model problem and that The- 



orem 4-8 applies to more general non-linearities. For example, the strong convex 
hull property applies to the non-linear problems mentioned in Remark \4-S\ if J is 
additionally Frechet differentiable, i.e., if the minimising problem \l.l\ can be equiv- 
alently formulated by a partial differential equation, the so-called Euler- Lagrange 



equation; compare with (1.2 I. 
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